Existence and stability analysis of finite O-tt-O Josephson junctions 
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We investigate analytically and numerically a Josephson junction on finite domain with two tt- 
discontinuity points characterized by a jump of tt in the phase difference of the junction, i.e. a 
O-TT-O Josephson junction. The system is described by a modified sine-Gordon equation. We show 
that there is an instability region in which semifluxons will be spontaneously generated. Using a 
Hamiltonian energy characterization, it is shown how the existence of static semifluxons depends 
on the length of the junction, the facet length, and the applied bias current. The critical eigenvalue 
of the semifluxons is discussed as well. Numerical simulations are presented, accompanying our 
analytical results. 

PACS numbers: 74.50.-|-r, 85.25.Cp, 74.20.Rp 
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I. INTRODUCTION 

Josephson junctions consist of two superconductors 
separated by a thin insulating barrier. When the in- 
sulation is thin enough, a current can flow across the 
barrier even when there is no potential difference. Cur- 
rent technological advances can manipulate the flow of 
the supercurrent such that its direction can depend on 
the spatial variable. This is possible, e.g., by imposing 
shifts to the Josephson phase. 

The idea of having a shift in the gauge phase 
of a Josephson junction was first proposed by Bu- 
laevskii et al,^ It was proposed that the presence 
of magnetic impurities may create a 7r-shift to the 
Josephson phase, which has been confirmed recentlyi^ 
Presently, one can also impose a tt phase-shift in a 
long Josephson junction using superconductors with 
unconventional pairing symmetry;^ Superconductor- 
Ferromagnet-Superconductor (SFS) 7r-junctions,^, 
Superconductor-Normal metal-Superconductor (SNS) 
junctionSfi or using a pair of current injectors^ All these 
findings have promising applications in information 
storage and information processing.^ 

This system, in which neighboring facets of a Joseph- 
son junction can be considered to have opposite signs of 
the critical current, present intriguing phenomena such as 
the intrinsic frustration of the Josephson phase over the 
junction and the spontaneous generation of a fractional 
magnetic fiux at the discontinuities i.e. the position of 
the jump in the Josephson phase i^iiiS 

In the present work, we consider the so-called O-tt-O 
Josephson junctions on a finite domain, modeled by a 
modified sine-Gordon equation with phase shift of = tt 
in some region and zero otherwise. An infinite domain 
O-TT-O Josephson junction was first studied by Kato and 
Imada,-^^ where they showed that there is a stability win- 
dow for the TT junction length in which the zero con- 
stant solution is stable. In the instability region, the 
ground state is non constant solution in space, which 
corresponds to antiferromagnetically ordered semiflux- 
ons. Later, using an existence analysis, it was shown that 
the presence of a minimum facet length of the tt junction. 



above which a non-trivial ground state emerges, also cor- 
responds to the minimum facet length needed to con- 
struct such solutions ji^ii^ The possibility of employing 
O-TT-O junctions for observations of macroscopic quantum 
tunneling was discussed in length by Goldobin et al.^^ In 
the presence of an applied bias current, a O-tt-O Joseph- 
son junction has a critical current above which one can 
thermally flip the order of the semifluxons^^ and another 
critical current above which the junction switches to the 
resistive state»i'^ Goldobin et a L^^d^ have also broadened 
the study of O-tt-0 junctions to 0-k-O junctions, where 
< K < TT (mod 27r). Here, we limit ourselves to discuss 
O-TT-O junctions only, but extend it to the case of a finite 
domain. This is of particular interest, especially from 
the physical point of view, as such junctions have been 
successfully fabricated recently,^^'^* which are certainly 
of finite length. 

The present paper is structured as follows. In Section 
II, we discuss the mathematical model that we use to 
describe the problem. We then show in Section III that 
when there is no bias current, the equation has two con- 
stant solutions. Due to the phase-shifts, there will be a 
region of facet lengths, in which both constant solutions 
arc unstable. In this instability region, a non-constant 
ground state will emerge from the constant solutions, i.e. 
a pair of semifluxons is the ground state of the system. A 
Hamitonian analysis is performed in Section IV to study 
the behavior of the non-trivial ground state, both with 
and without the presence of applied bias current. We 
then compare our analytical results with numerical com- 
putations in Section V. Finally, conclusions are presented 
in Section VI. 



II. MATHEMATICAL MODEL 

The dynamics of a finite Josephson junction with tt- 
discontinuity points is commonly described by the fol- 
lowing perturbed sine-Gordon equation 

(l)xx-(l>tt=sm{(l) + e{x))-^ + a(f>t, -L<x<L, (1) 

where a is a dimensionless positive damping coefficient 
related to quasi-particle tunneling across the junction, L 
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is the length of the junction, and 7 is the apphed bias cur- 
rent density normahzed to the junction critical current 
density Jc- Without loss of generality, in the following 
we set a — Q. 

Equation ([T]) is written after rescaling where the spa- 
tial variable x and time variable t are normalized to the 
Josephson penetration length Aj and the inverse plasma 
frequency ujp"^ , respectively. 

The function 6, representing the presence, or absence, 
of the additional 7r-phase shift, is given by 

0' L>\x\>a, 
^ ' \tt, \x\ < a, ^ ' 

where a is the length of the tt junction, i.e. the facet 
length. 

Equation ([1]) is subject to the continuity and boundary 
conditions 

(bi±a~) = 0(±a+), M±a-) = M±a+), M±L) = 0. 

(3) 

The governing equation ([T|), subject to the boundary 
conditions ([3]), can be derived from the Lagrangian 

C = J ^{\4>l -\<Pl-l + cos(0 + 9)+ -i<j,) dx. (4) 

As we mainly consider static semifluxons, the exis- 
tence of the solutions will be studied through the time- 
independent version of ([1]), namely 

(t>xx = sin(0 + 6*) - 7. (5) 

III. EXISTENCE AND STABILITY ANALYSIS 
OF CONSTANT SOLUTIONS 

It is clear that equation ([5]) admits two constant solu- 
tions (modulo 27r), namely 

(j) = arcsin7, tt — arcsin7, 

for a < < L, and 

4> — arcsin7 — tt, — arcsin7, 

for < < a. 

As solutions on the whole domain must satisfy the con- 
tinuity conditions Q, it can be concluded that constant 
solutions exist only when 7 = 0, i.e. (/) = and (jj — tt. 

Next, we will determine the linear stability of the con- 
stant solutions. For this purpose, we substitute the sta- 
bility ansatz 

(j) = 4> + ee^*V{x) (6) 

into Eq. ([1]). Neglecting higher order terms, one obtains 
the eigenvalue problem 

^ {E + co8{4, + 0))V, (7) 



where E = and V is also subject to the continuity and 
boundary conditions 

Vi±a-) = V{±a+), 14(±a-) = 14(±a+), 14(±L) = 0. 

(8) 

From the stability ansatz it is clear that (f) is stable 
if A has zero real parts. Therefore, a solution (j) is said to 
be linearly stable if < and unstable when E > 0. 

Due to the finite size of the domain, the eigenvalue 
problem ^ will give two sets of eigenvalues, i.e. in- 
finitely many eigenvalues that will constitute the contin- 
uous spectrum in the limit L 00 and a finite number 
in the discrete spectrum. For simplicity, in the following 
sections we denote the infinitely many eigenvalues by the 
'continuous' spectrum. 

A. Linear stability of 

First, we discuss the 'continuous' spectrum oi (p — 0. 
Looking for a bounded solution to ([7]) that satisfies the 
boundary conditions at a; = ±L, we obtain the solutions 

{A cos(Q;(a; + L)), ~L < x < —a, 
BCOS0X) + Csin(/3x), \x\ < a, (9) 
D cos{a:{x — L)), a < X < L, 

where a = V— 1 — E and (3 = %/! — E. 

Using the continuity conditions ([5]) , we obtain a system 
of four equations with four unknowns, given in a matrix 
form by 

/^\ 

Ml ^ = 0, 

\d) 

with a coefficient matrix Mi given in the appendix. 

To obtain a nontrivial V, we require det(Mi) = 0. An 
implicit plot of the equation, i.e. the 'continuous' spec- 
trum E(a,L), for L = 1 is shown in the top panel of 
Fig. [TJ Numerically, it is found that there is no unstable 
eigenvalue in the 'continuous' spectrum, i.e. E < for 
any a. As L increases, the distribution of E will become 
dense, as expected. 

Next, we find the discrete spectrum of the constant 
solution = 0, corresponding to bounded and decaying 
solutions of the eigenvalue problem ([7]). We obtain the 
solution 

{v4cosh(7(a; + i)), —L < x < —a, 
Bcos{Px)+Csin{Px), \x\ < a, (10) 
D cosh(7(a; — L)), a < x < L, 

where 7 — vT+~E and /3 is defined above. 

From the continuity conditions, again we find a system 
of four equations. As above, the eigenvalues are obtained 
by setting the determinant of the coefficient matrix M2, 
given in the appendix, to zero. An implicit plot of the 
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FIG. 1: Plot of the 'continuous' spectrum (top) and the dis- 
crete spectrum (bottom) of = as a function of the tt 
junction length for L — 1. 



eigenvalues as a function of a for L = 1 is shown in the 
bottom panel of Fig. [TJ 

From Fig. [Tl we observe that for a given L, there is 
a critical a above which E becomes positive, i.e., (f> = 
becomes unstable. In the following, we denote such a 
critical a by ac,o- For L — 1, ac,o ~ 0.46. As L increases, 
ac,o will asymptotically approach j, which is the critical 
facet length in the infinite domain calculated in Refs. 
[11II13I . The relation between Uco and L is implicitly given 
by the smallest positive root of 

cot(ac,o) tanh(i — ac,o) = —1, (11) 

which is obtained by considering the even mode of 
(fTU)) . i.e. setting E = C = and A^D = 
B cos{acfl) / cosh(L — ac,o)- For small L the root can be 
approximated by 

ac,o = ^-^L' + OiL'). (12) 

Plots of Qcfi as a function of L given implicitly by (|lip 
and its approximation (|12p are shown in Fig. [3) 
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FIG. 2: The same as Fig. [T] but for 4> = tt. 

B. Linear stability of tt 

Following the same steps as we did in the stability anal- 
ysis of = 0, the solution V to the eigenvalue problem 
([7]) that corresponds to the 'continuous' spectrum is given 
by 

{Acos{(3{x + L)), -L < X < -a, 
Bcos{ax) + C siii{ax), \x\ < a, (13) 
Dcos{l3{x - i)), a < X < L. 

One then finds that the spectrum is given by the zero 
of the determinant of the coefficient matrix M3, given in 
the appendix. 

A plot of the 'continuous' spectrum in the {E, a)-plane 
is shown in Fig. [21 from which it is clear that the contin- 
uous spectrum also only consists of stable eigenvalues. 

For the discrete spectrum of = tt in a finite do- 
main, bounded and decaying solutions V to the eigen- 
value problem ([7]) are given by 

{Acoa0{x + L)), -L < X < -a, 
Bcosh(7a;) -I- Csinh(7a;), \x\ < a, (14) 
D cos(/3(a; - L)), a < x < L. 

Due to the boundary conditions ([5]), again we obtain a 
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FIG. 3: Instability region of the constant solutions = and 
(j) = TT. Solid boundary curves are given by Eqs. pip and (|15|l . 
Dashed lines are analytical approximations given by (|12|) and 



system of four homogenous equations with a coefficient 
matrix M4, given in the appendix. 

The bottom panel of Figure [2] shows the plot of the 
zeros of det(M4) in the (iJ, a)-plane, for L = 1. We 
observe that for a close to zero, i?>0, i.e. = 7ris 
unstable. Yet, there is a critical value of a above which 
TT is stable. We denote this critical facet length by Qc-n, 
which for L = 1 is approximately 0.54. 

Again, considering the even state of pi)) . i.e. E — C = 
and A = D = i?cosh(ac,,r)/cos(L — ac,,r), one can show 
that ac,7r is related to L by the implicit equation 

coth(ac,7r) tan(L — ac,-n) — 1, (15) 
which, for L close to 0, can be approximated by 



(16) 



Combining Eqs. (fTTjl and (fTS]) . we plot in Fig. [3] the 
region in which both the stationary solutions of the sine- 
Gordon ^ are unstable. In the instability region, the 
ground state will be non constant in space. 



C. Symmetry 

Comparing Figs. [T] and [2 we observe that they are 
the same by reflection with respect to the line a = £/2, 
i.e. the stability of (/> = at the tt junction length a is 
the same as the stability of (/i = tt at the tt facet length 
(L—a). This symmetry occurs because, for the particular 
solutions, our Neumann boundary conditions at x = ±L 
^ can be replaced by periodic boundary conditions 



-L)=^{L), 0,(-L) =0,(L). 



(17) 



For the periodic boundary conditions, the governing 
equation ^ is symmetric by rotation, i.e. cyclic sym- 
metry, and 9 ^ 9 + TT. Using the symmetry, one can also 
conclude that 



for any L. 



IV. GROUND STATES IN THE INSTABILITY 
REGION 

In the following, we analyse perturbatively the ground 
states of the Josephson junction in the instability region. 
Our analysis, based on an Euler-Lagrange approxima- 
tion, is carried out for a facet length a close to one of the 
critical facet lengths ac,o and ac,-K- 



A. The case of < a — ac,o «C 1 
1. Existence analysis 

For a close to Oco, we approximate (/'(x) by 
cos(ac,o) 



[x] = B <^ cosh(L - Ucfi) 
cos(a;), \x\ < a 



cosh(a; -t- L), a < < L, 



(18) 

where B = B(t) and ac,o is given in (jlip . This expression 
of (j) is an exact solution to the linearization of ([5]) for a = 
flco and an arbitrary parameter B, i.e. we approximate 
the ground states by (/> = V{x) , where V{x) is the first 
even state to the eigenvalue problem ^ at E — and 
a = ac,o- 

Substituting the ansatz into the Lagrangian ([l]), 
writing a — L/2 — kL^ k < 1/24, and expanding about 
L = yields 



C = iL- L^/A)Bl - H, 
where the subscript represents a derivative and 
L'^B^ ( . ^2 



(19) 



H 



12 



24A:- 1 



^2L\L^(^-\k) + \-B-i 



-(6fc- 1) 



(20) 



The Euler-Lagrange equation from the Lagrangian 
(HH) is dt{dB,C) - deC = 0, giving 



Btt = 



1 



2i-L3/2 



(21) 



The time independent solution B — Bq of the Euler- 
Lagrange equation (pi]) is given by a cubic equation 
Hb = 0, or 



2L^Bo{B^{12k -2) + 3- 72k) 
9(L2 - 8) ■ 



(22) 
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For a general value of 7 ^ 0, we solve the cubic equation 
using Nickalls' method^ to obtain 

S^"' = 2Acos(e + 2(n- l)7r/3), n = l,2,3, (23) 
where 



A 



\Jk^, = arccos(-yAr//i)/3, 



yN = {2L ~ Ly4)-f, h = - 24A:)A. 

(n) 

When 7 = 0, the expressions for Bq are simplified to 



(1.2) 



/3(24fc- 



2(6fc- 1) 



(24) 



(12) 

The non-zero roots Bq ' represent a pair of antiferro- 
magnetically ordered semifluxons. 

If we study further, the three roots we find that 
they do not persist for all 7. If 7 is decreased (increased) 
away from zero, then there is a critical value of the bias 
curr ent at which B^^^ (S^^') colhdes with S^^^ in a saddle 
node bifurcation. Here, we denote this critical value of 
7 by 7c, 1- From our current approximation, 7c,ican be 
calculated from the condition y'j^ — h?^ which gives 



7c, 1 



2V2£^(24fc- 1)3/2 
9V6fc-l(L2-8) 



(25) 



2. Stability analysis 

To study the stability of the stationary solutions 

we easily check that when k < 1/24, H is locally mini- 

(1 2) 

mized by • To obtain the critical eigenvalue of the 

stable solutions, we write B = b'^^ + eB and substitute 
it into the Euler-Lagrange equation (|2ip to obtain 

1 



B{t)^ 



2i- £3/2 

The critical eigenvalue of B^^^ is then given by 

^=(2L-L3/2)^^^l5=<'' 
i.e. the square of the oscillation frequency of B{t). 



(26) 



(27) 



B. 



The case of < Oc 



a < 1 



To discuss the existence and stability of ground state 
solutions when a is close to ac,7r, we repeat the above cal- 
culations. We exploit the symmetry discussed in IIII CI 
that for the non-constant ground state, the Neumann 
boundary conditions ([3]) can be replaced by periodic 



boundary conditions (|17|) . Using this symmetry, we ob- 
tain that if 4'{x] a) is a ground state solution of the sine- 
Gordon equation with the tt facet length a, then 



{x; L — a) — TT — (f>{L ~ \x\; a). 



(28) 



Thus, the stability of the ground state in the limit 
< Uc^TT — a <C 1 can be deduced using this symmetry 
argument. 
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FIG. 4: (Top) Plot of <t>{0) and (t){±L) of the non-constant 
ground state obtained from numerical calculations (solid 
lines) as a function of the facet length a. Comparison with 
our analytical approximations (dashed lines) is also presented. 
(Bottom) The critical eigenvalue E of the solution depicted 
in the top panel. 



V. DISCUSSION 

To check our analytical results, we perform numeri- 
cal calculations and simulations. We numerically solve 
the time-independent governing equation ([5]), subject 
to boundary conditions ([3|) using a Newton-Raphson 
method, where we discretize the problem using central 
differences with a relatively fine spatial discretization. 
To numerically study the stability of a solution, we then 
discretize and solve the corresponding linear eigenvalue 
problem (cf. Eq. Q). 
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First, we study the existence and the stabihty of the 
non zero ground state in the absence of an apphed bias 
current. 

In the top panel of Fig. [H we plot (f>{x), which is rep- 
resented by 0(0) and 0(±L), of the non-constant ground 
states as functions of a for 7 = and L = 1. Due to 
the cyclic symmetry discussed in Section IIII CI we ob- 
serve that the curves are symmetric under rotation by 
TT radians. In the same figure, we also depict our ana- 
lytic approximation (j24p , where one can see a rather good 
agreement for a near to and Uc^tt- The approxima- 
tions deviate rapidly as a moves away from these critical 
lengths due to the fact that our junction length L is of 
order one. It is expected that for L close to zero, the 
approximations will give a rather good agreement. 

In the bottom panel of the same figure, we depict 
the critical eigenvalue of the non-zero ground states pre- 
sented at the top panel. At the critical facet lengths, the 
eigenvalues are certainly zero due to the change of stabil- 
ity with the constant solutions = 0, tt. Our analytical 
approximation (|27p is plotted in the same figure, from 
which we see that when the facet length a close to one of 
the critical values Ucfl and Uc^tt, the numerics are indeed 
well approximated by our analytical result. 

Next, we study the influence of an applied bias current 
to the existence and the stability of the non constant 
ground state. In the following, we particularly consider 
a = 0.495 and without loss of generality the 'positive' 
ground state, indicated by 0(0) > 0. The case of nega- 
tive 0(0) can be obtained simply by reflection due to the 
symmetry — > —(f> and 7 — s- —7. 

In the top panel of Fig. [5l we plot our numerical 0(0) 
as a function of the applied bias current 7 for L = 1 and 
a — 0.495. We use a path-following method starting from 
7 = 0. 

First, we decrease the applied bias current. As 7 is 
reduced, the value of 0(0) also decreases up to a certain 
value of bias current; the solution cannot be continued 
further, it terminates in a saddle node bifurcation. Using 
our path following method, the saddle-node bifurcation 
is indeed due to a collision with a non-constant solution 
bifurcating from = 0, as predicted by our analytical 
result. The value of 7 at which the bifurcation occurs is 
the aforementioned 7c, 1 (|25l) . Comparisons between the 
numerics and the analytical results of 7c, 1 are depicted 
at the middle panel of Fig. [51 

Besides decreasing 7, one can also increase it. As 7 
increases, the value of 0(0) also increases. As the bias 
current is increased further, a saddle-node bifurcation 
occurs. We denote this critical value of bias current by 
7c, 2- For a = 0.495, |7c,i| < 7c, 2- Using our path fol- 
lowing algorithm, we can follow the upper branch of the 
bifurcation and deduce that it corresponds to a collision 
between the non constant solution and = tt. Using the 
cyclic symmetry argument, we explain the bifurcation us- 
ing our analytical results Plotted in the top panel 
of Fig.[n]is our 7(^0) given by Eq. (P^ . properly shifted 
by TT, for the facet length {L~ a). The influence of the tt 
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FIG. 5: (Color online) The top panel depicts the existence 
diagram of the ground state. Plotted is 0(0) as a function 
of 7, obtained from numerical computations (solid lines) for 
a — 0.495 and L = 1. Shown in red is 0(0) as a function of a 
that corresponds unstable solutions. The upper and lower red 
branch corresponds to solutions tf) — tt, 0, respectively. The 
middle panel shows the critical bias currents 7c, i and 7c,2 as a 
function of a for L = 1. The bottom panel presents the critical 
eigenvalue of the non-constant ground state as a function of 
7 for a = 0.495 and L — 1. Analytical approximations are 
also presented in dashed and dash-dotted lines. 
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In the bottom panel of Fig. [5l we plot the critical eigen- 
value of the non constant ground state as a function of 7 
for a = 0.495 and L = 1. Interestingly, we observe that 
the lowest eigenvalue is attained at a non-zero bias cur- 
rent. This indicates that a non constant ground state can 
be made to be 'more stable' by applying a bias current. 
On the same figure, we also plot our approximations (|27p . 
which qualitatively agree with the numerical results. 

Studying further the saddle-node bifurcation between 
the non constant ground state and (/> = tt in Fig. [SJ we 
observe that it is not the typical collision that leads to 
the definition of 7c. 2 for any L. When L is relatively 
large, we find that the upper branch does not necessarily 
correspond to a constant solution. In Fig. [6l we consider 
another case for L = 10 and a = 3. 

Starting on the middle branch from 7 = and 0(0) « 
2.5, we then increase the bias current. At the critical bias 
current 7^2, i-e. 7 ~ 0.6, we have a saddle-node bifurca- 
tion. Using our path following code, we follow the branch 
beyond the bifurcation point, from which we obtain that 
the branch does not correspond to a constant solution. 
In the middle panel of Fig. [51 we plot the corresponding 
solutions for some values of 7. Considering the profile 
(/>(x) at 7 = 0, we could conclude that it corresponds to a 
pair of semifluxons and one fluxon on each side. The pro- 
file is similar to the so-called semifluxon type 3, defined 
in Ref. Is' for an infinitely long O-tt Josephson junction. 
From our numerical computations (not shown here), the 
two branches seem to be distinguished by the ability of 
the junction of length L to support an additional fluxon 
on both sides. The first critical value 7c, 1 corresponds to 
the collision between the non constant ground state and 

= 0. The bottom panel presents 7c. 1 and 7c, 2 as a func- 
tion of a for L = 10. An approximate expression for the 
critical currents is presented in dashed lines given 
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FIG. 6: (Color online) Top panel is the same as the top panel 
of Fig.O but for L = 10 and a = 3. Middle panel shows some 
of the corresponding solutions of the top branch for different 
values of 7. Bottom panel shows 7c, 1 and 7c, 2 as a function 
of a. Solid and dashed lines are numerical calculations and 
approximations 



facet length on the existence diagram is indicated by the 
arrows, i.e. as a increases (decreases) towards ac,-n- the 
two lobes move according to the arrows (and vice versa). 
In the middle panel of Fig. [5l we plot the second critical 
bias current, 7c, 2, as a function of the tt facet length a. 
Again, due to the cyclic symmetry, 7c. 2 can be obtained 
from 7c. 1 by rotating the curve by tt radians. 



where ac,o ~ 7r/4 and a^.-K ~ {L — n/A). 

One may ask about using O-tt-O Josephson junctions to 
observe macroscopic quantum tunneling. To answer the 
question, the reader is addressed to Ref. which con- 
siders quantum tunneling of a semifluxon in a finite O-vr 
junction, where it was concluded that finite O-tt junctions 
do not provide a good playground to build a qubiti^ Us- 
ing the similarity between the currently used Neumann 
boundary conditions and periodic boundary conditions 
(see Section IIII Cp , one may conclude that the system 
considered here is also not a promising one in which to 
observe macroscopic quantum tunneling, as finite O-tt-O 
junctions considered here are made of two finite O-tt junc- 
tions. 
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VI. CONCLUSION 

We have investigated analytically and numerically O-tt- 
Josephson junctions on a finite domain. Wc have shown 
that there is an instability region for constant solutions 
in which semifluxons are spontaneously generated. Us- 
ing an Euler-Lagrange approximation, it has been shown 
that the existence of static semifluxons depends on the 
length of the junction, the facet length, and the applied 
bias current. In addition the critical eigenvalue of the 
semifluxons has been discussed. Numerical simulations 
have been presented, accompanying our analytical re- 
sults. 

For future investigations, the two dimensional version 



of Josephson junctions with phase-shifts, = tt in some 
areas and 6 = elsewhere, will be considered. The effects 
of boundary conditions on the stability of non constant 
ground states in such a system will certainly be of inter- 
est. These are works in progress and will be reported in 
future publications. 



APPENDIX A: COEFFICIENT MATRICES 

The coefficient matrices Af„. n = 1,2,3,4, used to 
derive the 'continuous' and the discrete spectrum of 
= 0, TT are given by 
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Q;sin(Q!(L — a)) 





— cos (/3a) 
-0sm{$a) 
cos(/3a) 
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—$cos{$a) 

sin(/3a) 
/3cos(/3a) 






— cos(Q!(a — L)) 
Q;sin(d(a — L)) 



(Al) 



Mo 



cosh(7(L — a)) — cos(/3a) sin(/3a) 

7sinh7(L — a) —$sm{$a) — /3cos(/3a) 

cos(/3a) sin(;9a) cosh(7(a — L)) 

—$sm{(}a) /3sin(/3a) 7sinh(7(a — L)) 



(A2) 
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— dsin(da) $sm{(3a) — /3sin(/3(o — L)) 



(A3) 
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